arXiv:1501.03565vl [astro-ph.HE] 15 Jan 2015 


Mon. Not. R. Astron. Soc. 000,[T]|3() Printed 16 Januaiy 2015 (MN MgX style file v2.2) 


A disc corona-jet model for the radio/X-ray correlation in black hole 
X-ray binaries 


Erlin Qiao and B.F. Liu ^ 

'National Astronomical Obsen’atories, Chinese Academy of Sciences, Beijing 100012, China 


16 January 2015 


ABSTRACT 

The observed tight radio/X-ray correlation in the low spectral state of some black hole X- 
ray binaries implies the strong coupling of the accretion and jet. The correlation of Lr oc 
^ was well explained by the coupling of a radiatively inefficient accretion flow and 
a jet. Recently, however, a growing number of sources show more complicated radio/X-ray 
correlations, e.g., Lr oc for LxIL^a <: 10 “^, which is suggested to be explained by 
the coupling of a radiatively efficient accretion flow and a jet. In this work, we interpret the 
deviation from the initial radio/X-ray correlation for Lx/TEdd ^ 10“^ with a detailed disc 
corona-jet model. In this model, the disc and corona are radiatively and dynamically coupled. 
Assuming a fraction of the matter in the accretion flow, rj = Mjet/M, is ejected to form the 
jet, we can calculate the emergent spectrum of the disc corona-jet system. We calculate Lr 
and Lx at different M, adjusting 77 to fit the observed radio/X-ray correlation of the black hole 
X-ray transient H1743-322 for Lx/Ledd > 10“^. It is found that always the X-ray emission is 
dominated by the disc corona and the radio emission is dominated by the jet. We noted that 
the value of 77 for the deviated radio/X-ray correlation for Lx/Ledd > 10“^, is systematically 
less than that of the case for Lx/Lndd < 10“^, which is consistent with the general idea that 
the jet is often relatively suppressed at the high luminosity phase in black hole X-ray binaries. 

Key words: accretion, accretion discs — black hole physics — X-rays; individual: H1743- 
322 — X-rays: binaries — radio continuum: stars 


1 INTRODUCTION 

A black hole X-ray binary (BHB) is a gravitationally bound sys¬ 
tem composed of a black hole and a normal star. BHBs are lumi¬ 
nous in X-rays, which is believed to be resulted by accreting the 
matter of the normal star onto the black hole. According to the X- 
ray spectral features and the timing properties, two typical spectral 
states were identified in BHBs, i.e., the high/soft spectral state and 
the low/hard spectral state (for reviews, see Tanaka & Lewin 1995; 
Remillard & McClintock 2006). In the high/soft spectral state, the 
spectrum is dominated by a peak emission around IkeV, which is 
believed to be produced by a cool disc extending down to the in¬ 
nermost stable circular orbit (ISCO) of a black hole (Shakura & 
Sunyaev 1973; Mitsuda et al. 1984; Belloni et al. 2000). The spec¬ 
tral features of the low/hard spectral state are complicated, which 
are generally thought to be produced by an inner hot accretion flow 
and an outer truncated cool disc (Rees et al. 1982; Esin et al. 1997, 
2001; McClintock et al. 2001; Yuan et al. 2005; Narayan & Mc¬ 
Clintock 2008). Recently, some observations indicate that a cool 
disc may also exist in the region very close to ISCO in the low/hard 
spectral state (Miller et al. 2006a, b; Tomsick et al. 2008; Reis et 
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al. 2010). Meanwhile, the low/hard spectral state is often associ¬ 
ated with the production of collimated, relativistic jets, which are 
quenched in the high/soft spectral state (Fender et al. 2004). 

A correlation between the radio luminosity and X-ray lu¬ 
minosity was found in the low/hard spectral state of BHBs, i.e., 
Lr oc L* , with ~ 0.5 - 0.7 (Corbel et al. 2003, 2008, 2013; Gallo 
et al. 2003). The existing of the radio/X-ray correlation for BHBs 
is mainly from the observations of two sources, i.e., GX 339-4 and 
V404 Cyg (Corbel et al. 2013). Yuan & Cui (2005) interpreted 
this radio/X-ray correlation, i.e., Lr oc for Lx - 
to Lx - lO^^Ledd (with Lndd = 1.26 x 10 ^*M/Mo erg s“' ) within 
the framework of a radiatively inefficient accretion flow (RIAF)-jet 
model, in which a fraction of the matter, 77 , in the accretion flow 
is assumed to be ejected to form a jet (with 77 = Mjet/M). In this 
model, the radio emission was dominated by the self-absorbed syn¬ 
chrotron emission of a steady, collimated compact jet, and the X- 
ray emission was dominated by a RIAF via thermal Comptoniza- 
tion process. In the RIAF-jet model, if a constant 77 is assumed 
for different M, the radio luminosity from the jet can be scaled as 
Lr oc t]M^ with ^ ~ I (e.g., Heinz & Sunyaev 2003), and because 
of the nature of the advection in RIAF, the X-ray luminosity can 
be scaled as Lx 00 M'? with q ~ 2. Then, the predicted radio/X-ray 
correlation is Lr oc L^^'' with^/q ~ 0.5, which is roughly consistent 
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with observations. The RIAF-jet model can not only interpret the 
radio/X-ray correlation, but also can explain the broadband spectral 
energy distribution (SED) and most of the complex timing features 
of some black hole X-ray transients, e.g., XTE Jill 8+480 (Yuan et 
al. 2005; Malzac et al. 2004). When BHBs enter the quiescent state, 
i.e., Lx Yuan & Cui (2005) predicted a steeper radio/X- 

ray correlation of Lr oc within the framework of the RIAF-jet 
model, in which the radio emission is dominated by the jet, mean¬ 
while the X-ray emission is also dominated by the jet rather than the 
RIAF. Later, by fitting the SEDs of the black hole X-ray transient 
GRO J1655-40 and V404 Cyg in the quiescent state, it is found that 
both the radio emission and the X-ray emission can indeed be ex¬ 
plained by the jet (Pszota et al. 2008; Xie et al. 2014). However, 
due to the very faint emission in the quiescent state, the presence of 
the radio/X-ray correlation of Lr oc for Lx < lO^^’LEdd is still 
controversial (Gallo et al. 2006). 

We need to keep in mind the complexities of the observed 
radio/X-ray correlation. For example, by critically examining the 
radio/X-ray correlation in a sample of microquasars, Xue & Cui 
(2007) found that the correlation varied significantly among indi¬ 
vidual sources, not only in terms of the shape but also of the degree 
of the correlation. Recently, a statistical analysis of the data led to 
the new claim of dual tracks, with some referring as the ‘univer¬ 
sal’ track and the other as being ‘outlier’ track for LxIL^a <: 10 “^ 
(Gallo et al. 2012; Jonker et al. 2012). We should note that the ob¬ 
served radio/X-ray correlation is quite complex, and the claim of 
the simple two tracks is also fairly dubious (Coriat et al. 2011). 
Anyhow, so far, a growing number of sources have been discov¬ 
ered with a steeper radio/X-ray correlation, namely, Lr oc Lj/' ‘*, 
during Lx ^ lO'^LRd,) [e.g., IGR J17497-2821 (Rodriguez et al. 
2007), XTE J1650-500 (Corbel et al. 2004), Swift J1753.5-0127 
(Soleri et al. 2010)]. The study of the radio/X-ray correlation of 
Lr oc * for Lx <: lO^^LR^d is the purpose of this work. By 

collecting the archive quasi-simultaneous radio and X-ray data be¬ 
tween 2003 to 2010, Coriat et al. (2011) comprehensively studied 
the relationship between the radio luminosity and X-ray luminosity 
of the black hole X-ray transient H1743-322. It is found that during 
a high luminosity phase, corresponding to Lx/Ledd ~ 10 “^ to 10 “', 
H1743-322 follows the radio/X-ray correlation with a steeper slope 
ofb~ 1.4. While during a low luminosity phase, corresponding to 
Lx/Lsdd £ 10“^, H1743-322 follows the radio/X-ray correlation 
with a slope of b ~ 0.6. For LxIL^a between 10“^ to 10“^, it is 
probably corresponding to a transition region between the two cor¬ 
relations. 

The change of the radio/X-ray correlation from the low lu¬ 
minosity phase to the high luminosity phase in H1743-322 may 
imply either the change of the properties of the accretion flow or 
the change of the different coupling in the accretion flow and jet, 
e.g., the change of the dependence of the fraction of the matter rj 
ejected to form the jet on the mass accretion rate M. The sources 
with the radio/X-ray correlation of Lr oc Lj/'-'* could be considered 
as a X-ray loud hypothesis, i.e., for a given radio luminosity, the 
simultaneous X-ray luminosity of the track with b ~ 1.4 is higher 
than that of the track with b ~ 0.5 - 0.7, or in turn could be con¬ 
sidered as a radio quiet hypothesis, i.e., for a given X-ray luminos¬ 
ity, the simultaneous radio luminosity of the track with b ~ 1.4 is 
lower than that of the track with b ~ 0.5 - 0.7. In both hypotheses, 
the radio emission is dominated by the self-absorbed synchrotron 
emission of the steady, collimated compact jet, and the X-ray emis¬ 
sion is dominated by the accretion flow. The main difference of fhe 
two hypotheses is the different dependence of 77 on M. In the X- 
ray loud hypothesis, if a constant 77 is assumed for different M, for 


a given radio luminosity, the higher X-ray luminosity of the track 
with b ~ 1.4 is probably resulted by the transition of the accretion 
flow from a RIAF fo a radiatively efficienfly accretion flow, e.g, a 
disc corona system by Haardt & Maraschi (1991, 1993), Di Matteo 
el al. (1999), Liu el al. (2002a, 2003), Merloni & Fabian (2002), 
Cao (2009) and Huang el al. (2014), or a luminous hot accretion 
flow by Yuan (2001), Xie & Yuan (2012) and Ma (2012). Meyer- 
Hofmeister & Meyer (2014) proposed that recondensation of gas 
from the corona into an inner disc can provide additional soft pho¬ 
tons for Comptonization, which leads to a higher X-ray luminosity 
compared to the unchanged radio emission. In the radio quiet hy¬ 
pothesis, a varied 77 is assumed for different M. As suggested by 
Coriat et al. (2011), if there is a liner dependence of 77 on M, i.e., 
Tj ac M, the predicted radio/X-ray correlation should be Lr oc . 
Because for RIAF q is ~ 2, the predicted slope of the radio/X-ray 
correlation is also roughly close to the track with b ~ 1.4. How¬ 
ever, theoretically the dependence of 77 on M is unclear, (e.g., Pe’er 
& Casella 2009), further studies are still needed to put constraints 
on the relation between 77 and M. 

Observationally, there is evidence of a coupling of the hot 
plasma and the jet in both BHBs and active galactic nuclei (AGNs), 
i.e., the coupling of the disc corona and jet at high mass accretion 
rates and the coupling of the RIAF and jet at low mass accretion 
rates (Wu et al. 2013; Zdziarski et al. 2011; Merloni et al. 2003; Fal- 
cke et al. 2004). Meanwhile, theoretically, for M > (with a 

the viscosity parameter, M^i = 1.39 x 10**M/Mo g s“'), the accre¬ 
tion flow will transit from a RIAF to a disc corona system, which 
has been comprehensively studied by many authors, e.g., Meyer et 
al. (2000a, b), Liu et al. (2002b) and Qiao & Liu (2009, 2010, 2013) 
within the framework of the disc evaporation model, or by Narayan 
& Yi (1995b), Abramowicz et al. (1995) and Mahadevan (1997) 
within the framework of the RIAF solution. 

Consequently, in this paper, we propose a disc corona-jet 
model to explain the radio/X-ray correlation of Lr oc LJ* '' during 
the high luminosity phase > 10 “^, in which a fraction of 

the matter, 77 , in the corona is assumed to be ejected to form the jet. 
So far, the theoretical understanding of the jet formation is poor. 
Specifically, it is difficult to put constraints on the dependence of 
77 on Af in our model, so we set 77 as an independent parameter on 
M to fit the observations. As an example, by fitting the observed 
radio/X-ray correlation during the high luminosity phase of black 
hole X-ray transient H1743-322 for Lx/Ledd > 10“^, we found that 
77 is weakly dependent on M, and the mean fitting result of 77 is 
~ 0.57%. The derived relatively smaller fitting value of 77 supports 
the general idea that the jet is often suppressed at the high luminos¬ 
ity phase of BHBs. The disc corona-jet model is briefly described 
in section|2] The numerical results are presented in Section[3 Some 
comparisons with observations are shown in Section 4. Discussions 
are in Section 5, and the conclusions are in section 6 . 


2 THE MODEL 
2.1 Accretion flows 

The accretion flows adopted here is a geometrically thin disc en¬ 
closed by a geometrically thick, hot corona around a central black 
hole (Liu et al. 2002a, 2003). The disc is a standard Shakura & Sun- 
yaev (1973) disc, which is tightly coupled with the plane-parallel 
corona. Magnetic fields are assumed to be generated by dynamo 
action. As a result of Parker instability, magnetic flux loops contin¬ 
uously emerge from the disc to the corona and reconnect with other 
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loops. In this way, the accretion energy taken from the thin disc is 
released in the corona as thermal energy and eventually emitted 
away mostly in X-ray band via inverse Compton scattering. The 
density of the corona is determined by an energy balance between 
the downward heat conduction and mass evaporation in the chro¬ 
mospheric layer. A detailed description of the model can be found 
in Liu et al. (2002a, 2003). The equations describing these pro¬ 
cesses in the corona are listed as follows. 
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where T is the coronal temperature, n is the coronal number den¬ 
sity, B is the strength of magnetic field, U^d is the energy density 

- /z 


of the soft photon field for Compton scattering, Va = y ^ the 
Alfven speed, £ is the length of magnetic loop, r* is the effective 
optical depth, defined as t* = A^ncril with ~ 1. Other con¬ 
stants have their standard meanings. From equations (1) and (2), 
the temperature T and density n in the corona can be determined 
for a given magnetic field B and radiation field (/rad> and then the 
radiative spectrum can be calculated out. 

For a coupled disc and corona system, the magnetic field is 
derived from an equipartition of gas pressure and magnetic pres¬ 
sure in the disc, i.e., p = = 1- The soft photons (I/rad) 

are composed of intrinsic disc radiation and the reprocessed 
radiation of backward Compton emission of the corona With 
energy transporting to the corona by the magnetic field, the disc is 
expected to be gas pressure dominant. Thus, the magnetic field and 
soft photon field can be expressed as functions of black hole mass, 
accretion rate, and distance. 
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where m, mo.i, oio.i, P\, rio are the mass of black hole, the ac¬ 
cretion rate, the viscosity parameter, the equipartition factor, and 
the distance respectively scaled by Mq, 0 . IMndd, 0.1, 1 and lOf^s- 
0 = 1 - and R, - 3Rs is adopted as the ISCO of a non¬ 

rotating black hole {Rs= 2.95 x 10^ M/Mg cm). A„ is a factor in¬ 
troduced for the evaluation of the seed field in Haardt & Maraschi 
(1991, 1993), which is around 1 in order of magnitude. Here the 
energy fraction dissipated in the corona, /, is not a free parameter 
but can be expressed as. 
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Combing equations Q, 0, 0, 0 and 0, we get a solution 
of the corona above a gas pressure-dominated disc in the case of 
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Given the values of the input parameters m, m, a, £, and the 
initial values of Ar and d,,, we can solve equation 0 for /. Then T 
and n are solved from equations 0 and 0. The effective tempera¬ 
ture of the soft photos Tr 

_ 3 CMM ^ l - f } c,,re 

CIO) 

C T ' ' 

~ 4^rad 

can also be calculated by combing equations 0 and 0 (with 
albedo always being assumed to be zero in our calculations). With 
T, n and Tr, the spectra of the disc-corona system can be calculated 
by Monte Carlo simulation. 

In our detailed calculations, we also take into account the in¬ 
fluences of Ar and A„ and the length of magnetic loops. Since the 
corona temperature T, density n and energy fraction / do not sen¬ 
sitively depend on Ar and A„ (see equations 0-0), and the value 
of Ar and /i„ should be ~ 1 in order of magnitude, the corona spec¬ 
tra can not be significantly affected by the chosen value of Ar and 
A„. Nevertheless, we repeat the Monte Carlo simulation by adjust¬ 
ing the value of Ar and d,,. We find a set of reasonable value for 
Ar and A,, until the downward-scatted luminosity is equal to the 
seed luminosity and the upward-scatted luminosity is equal to the 
released gravitational energy. In this way we find a self-consistent 
solution of the disc-corona system and get the corresponding emer¬ 
gent spectrum. We test the effect of the length of magnetic loop £ 
on the emergent spectrum, it is found that the emergent spectrum 
is very weakly dependent on £. So in our model, £ is not a free pa¬ 
rameter, and we set £ = lOf^s throughout the calculations. For the 
detailed description of the Monte Carlo simulation, one can also 
refer to Liu et al. (2003). 


2.2 Coupled disc corona-jet model 

The calculation of jet emission is based on the internal shock sce¬ 
nario as described in Yuan et al. (2005). In the disc corona-jet 
model, a small fraction of the matter in the accretion flow, i.e., 
Mjet = r]M, is assumed to be ejected to form the jet. In this work, 
a conical geometry is considered for the jet, and the half opening 
angle of the jet is fixed al ip = 0.1. Changing the value of <p will 
lead to a change of the density of the jet, however this effect can be 
absorbed by the mass accretion rate. The bulk Lorentz factor of the 
jet is fixed at Tjet — 1-2, which is well consistent with the typically 
observed value of the jet in the low/hard spectral state of BHBs (e.g. 
Fender 2006). Within the jet, internal shock is produced due to the 
collision of the shells with different velocities. The internal shocks 
will accelerate a small fraction of the electrons to be a power-law 
energy distribution with index p. We fix p = 2.1 in the calculations 
as suggested by Zhang et al. (2010) for fitting the SEDs of three 
black hole X-ray transients J1753.5-0127, GRO J1655-40 and XTE 
J1720-318. The two parameters, and eb, describing the ratio of 
the energy of the accelerated electrons and the amplified magnetic 
field to the shock energy in the shock front are fixed at = 0.04 
and fB = 0.02 respectively, which are consistent with the obser¬ 
vations of gamma-ray bursts (GRB) afterglows (e.g., Panaitescu & 
Kumar 2001, 2002). Because the Compton scattering optical depth 
of the jet is small, only synchrotron emission is considered in the 
calculation (Markoff et al. 2001). 


3 NUMERICAL RESULT 

We calculate the emergent spectra of the disc corona-jet model 
around a stellar-mass black hole with mass M when the parameters 
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including M, a, and 77 are specified. In this paper, we fix the black 
hole mass at M = IOMq, assuming a typical viscosity parameter 
a = 0.3, and a constant 77 for different M. 

Given 77 = 0.2%, the emergent spectra are plotted in the left 
panel of Figure [T] for different mass accretion rates, i.e., from the 
bottom up M = 0.02,0.05,0.1,0.3 and 0.5 Medd- The solid lines 
are the total emergent spectra, and the dotted lines are the emergent 
spectra from the jets. From the left panel of Figure|T] we can clearly 
see that X-ray emission is dominated by the disc-corona system and 
the radio emission is dominated by the jet for the mass accretion 
rates from M = 0.02 - 0.5 Meaa- Based on the emergent spectra, 
we calculate the radio luminosity Lg.soHz and the X-ray luminosity 
T 2 -iokeV for different M. The best-fitting linear regression for the 
correlation between Tg.soHz/TEdd and M can be expressed as. 


Tg.SGHz I 

Lem ) 


= A(77) 
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L] 




Mem ) 


( 11 ) 


where A( 77 )|,=o, 2 % = lO^*'^^ and ^( 77 )|,,=o. 2 % = l-42. The best-fitting 
linear regression for the correlation between T 2 -iokev/TEdd and M 
can be expressed as, 


T2-10keV \ _ yq -1.23{ 

Lem ) l^Eddi ’ 
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where q = 1.06, which is roughly consistent with a radiatively ef¬ 
ficient accretion flow with Lx oc M (e.g. Haardt & Maraschi 1991, 
1993). Combing equations [TT]and[T2 we can derive that, 
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where ^( 77 )/^|,=o. 2 % ~ 1-35 and A( 77 ) 10 ‘'^^f'''*^'’|,=o. 2 % = 10“®*. The 
derived slope of radio/X-ray correlation from the disc corona-jet 
model is steeper than that of the RIAF-jet model, i.e., for a fixed 
X-ray luminosity, the increase of the radio luminosity of the disc- 
corona system is more quick than that of the RIAF. One of the 
reasons is that, for a given radio luminosity, due to the nature of the 
high radiative efficiency, the simultaneous X-ray luminosity pre¬ 
dicted by the disc-corona system is intrinsically higher than that 
of the RIAF. Moreover, because the X-ray luminosity of RIAF is 
roughly Lx or M^, and the X-ray luminosity of the disc-corona sys¬ 
tem is roughly Lx oc M* ®®, for a fixed radio luminosity, the increase 
of the X-ray luminosity of the RIAF is more quick than that of the 
disc-corona system, which in turn is equivalent to the case, i.e., for 
a fixed X-ray luminosity, the increase of the radio luminosity of the 
disc-corona system is more quick than that of the RIAF. 

In order to check the effects of 77 on the radio/X-ray correla¬ 
tion, we assume another constant 77 for different M, i.e., 77 = 0.5%, 
to calculate the emergent spectra for comparisons with that of 
77 = 0.2%. The emergent spectra are shown in the right panel of 
Figure Ul The best-fitting linear regression for the correlation be¬ 
tween Ls^GHz/LFdd and M can be expressed as, Lq^cmvILem = 
A(77)(M/MEdd)«''' with A( 77 )|,=o. 5 % = IQ-’’" and = 

1.40, and the best-fitting linear regression for the correlation be¬ 
tween Ld-iokev/Tfidd and M can be expressed as, L 2 _iokev/TEdd = 
{M!MemT with q = 1.06. Then, the predicted radio/X- 
ray correlation is Lg.SGHz/LEdd = A(77)10‘'2®«'')/'?(L2-10kev/LEdd)^‘''>/'' 
With A(77)10‘«fW/^|,.o.5% = 10-®* and ^(77)/g|,=o.5% ~ 1-32. 
Roughly, it is clearly shown that a systematic increase of the frac¬ 
tion of the matter ejected to the jet leads the radio flux to a system¬ 
atic increase, however does not change the slope of the radio/X-ray 
correlation much. 

In the above calculations, we assume that 77 is a constant for 


different M to produce the radio/X-ray correlation. However, actu¬ 
ally the dependence of 77 on Af is unclear, so 77 will be set as an 
independent parameter on M to fit the detailed observations in the 
next section. 


4 APPLICATION TO BLACK HOLE X-RAY TRANSIENT 
HI743-322 

H1743-322 is a X-ray transient discovered in 1977 (Kaluzienski & 
Holt 1977; Doxsey et al. 1977), then detected again by the Interna¬ 
tional Gamma-ray Astrophysics Laboratory (INTEGRAL) during 
the burst in 2003. H1743-322 is well studied in both radio and X- 
ray band. By comparing both the spectral features and the timing 
properties of H1743-322 to the well-studied black hole X-ray tran¬ 
sient XTE J1550-564, McClintock et al. (2009) argued that H1743- 
322 is a accreting black hole. The black hole mass of H1743-322 
is inferred to be IOMq (e.g., Steiner et al. 2012). Coriat et al. 
(2011) analyzed all the archive data of H1743-322 observed by 
RXTE between 2003 January 1 and 2010 Februay 13, meanwhile 
the authors conducted quasi-simultaneous (Af < Id) radio observa¬ 
tions through Telescope Compact Array (ATCA), Compact Array 
Broadband Backend (CABB), and collected the observational data 
of VLA in the literature. 

Based on the data of H1743-322 from Coriat et al. (2011), 
we plot the relation between Ls-gkev/LEdd and Lg.sGHz/LEdd with 
the red sign ‘•’in the right panel of Figure |2] For L 3 _ 9 i;ev/LEdd > 
10 “*, the best-fitting linear regression for the correlation between 
Ls-Qkev/LEdd and Lg.sGHz/LEdd is as follows, 

jQ-S.gf T3-9keV 
( LEdd 

This is plotted in the right panel of Figure |2] with the dotted line. 
Fixing M = IOMq, a = 0.3, we calculate the radio luminosity 
Lg.sGHz and X-ray luminosity L 3 _ 9 itev for different M, adjusting 77 
to fit equation [T4] It is found that 77 is very weakly dependent on 
M, i.e., for M = 0.02 MEdd^ 'I = 0.54%; for M = 0.05 MEdd^ 7 = 
0.52%; for M = 0.1 Mem,, 77 = 0.57%; for M = 0.3 Mem, J] = 
0.62%; for M = 0.5 MEdd^ U = 0.62%. The mean fitting value of 
77 is ~ 0.57%. The thick solid line in the right panel of Figure |2] 
is the model line. The corresponding emergent spectra with mass 
accretion rates are plotted in the left panel of Figure[2] Yuan & Cui 
(2005) fitted the radio/X-ray correlation of Lr oc LjJ® * for Lx - 
10“®LEdd to Lx - 10“*LEdd within the framework of the RIAF-jet 
model. They found that the fitting result was very sensitive to a 
parameter 6, which denotes the fraction of the viscosity dissipated 
energy directly heating the electrons in the RIAF. For 6 = 0.5, 
they found that 77 was highly dependent on M, decreasing from 77 ~ 
10% to 77 ~ 1% with mass accretion rate from M ~ 10-*MEdd to 
~ 10-*MEdd- For a smaller value of d = 0.01, they found that 77 
was nearly a constant with M, i.e., 77 ~ 1%. Indeed, we have tested 
the value of 77 for H1743-322 during the phase LxILem £ 10“* 
within the framework of RIAF-jet model. We found 77 10%. It 

is clear that the fitting value of 77 from the RIAF-jet model for the 
low luminosity phase is systematically higher than that of the fitting 
value of 77 from the disc corona-jet model for the high luminosity 
phase, which is consistent with the general idea that the jet is often 
relatively suppressed at the high luminosity phase in BHBs (Fender 
2004). 

In the present jet model, besides the mass rate in the jet Mjet, 
there are still three parameters, which can affect the emission of the 
jet, i.e., the power-law energy distribution index of the accelerated 
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Figure 1. 


Emergent spectra of the disc corona-jet model around a black hole with M = IOMq assuming a = 0.3. Left panel, r] = 0.2% is assumed. From the bottom up, 
the solid lines are the combined emergent spectra of the disc corona and jet for M = 0.02,0.05,0.1,0.3 and 0.5 Mem respectively. The dotted lines are the 
emergent spectra from the jets. Right panel: r] = 0.5% is adopted and the meaning of the line styles are the same as in the left panel. 


electrons p, the ratio of the energy of the accelerated electrons and 
the amplified magnetic field to the shock energy, and . In this 
work, p = 2.1, fe = 0.04 and = 0.02 are fixed respectively. 
By modeling the broadband emission of the afterglow of eight 
GRBs, Panaitescu & Kumar (2001) derived that p = 1.87 ± 0.51, 
Ce = 0.062 ± 0.045 and logen = -2.4 ± 1.2. Yuan et al. (2005) 
fitted the multiwavelength observations of XTE 11118+480 with 
RIAF-jet model, in which the best fitting results are p = 2.24, 
ec = 0.06 and es = 0.02 respectively. Based on RIAF-jet model, 
Zhang et al. (2010) fitted the simultaneous multiwavelength obser¬ 
vations of three black hole X-ray transients 11753.5-0127, GRO 
11655-40, and XTF 11720-318. It is found that the best fitting re¬ 
sults in 11753.5-0127 are p = 2.1, = 0.04 and eg - 0.02, in GRO 

11655-40 are p = 2.1, ee = 0.06 and eg = 0.02, and in XTF 11720- 
318 are p = 2.1, e,. = 0.06 and £3 = 0.08. It has been tested that the 
different values of p has very minor effects on the radio emission, 
however has significantly effects on the X-ray emission (Ozel et al. 
2000; Figure 4 in Yuan et al. 2003). Due to the X-ray emission is al¬ 
ways dominated by the emission of disc and corona, the change of p 
in an observationally reasonable range will not change our results. 
Throughout the calculation p = 2.1 is fixed. We test the effects of 
£5 and eg on the jet emission. In the left panel of Figure fixing 
M = IOMq, a = 0.3, M = 0.5, 77 = 0.5%, p = 2.1 and £3 = 0.02, 
we plot the emergent spectra of jet for e^ = 0.01,0.04,0.08,0.1 re¬ 
spectively. It is found that, e^ is sensitive to the X-ray emission, and 
insensitive to the radio emission. However, in the observationally 
inferred range of e^ ~ 0.01 - 0.1, it is clear that the X-ray emis¬ 
sion is still dominated by the disc and corona instead of the jet, so 
the change of £e will not change our results (see the right panel of 
Figure[T]for comparison). Throughout the calculation e^ = 0.04 is 
fixed. In the right panel of Figure[3 fixing M = IOMq, a = 0.3, 
M = 0.5, r] = 0.5%, p = 2.1 and £e = 0.04, we plot the emergent 
spectra of jet for £3 = 0.01,0.02,0.05,0.1 respectively. It is found 
that, increasing the value of £3 will increase the radio emission, and 
nearly does not change the X-ray emission. As discussed in Section 
3.1, an increase of Mjet will also increase the radio emission. So de¬ 
tailed spectral fitting are needed to disentangle the effects of £3 and 


Mjet on the jet emission, which is beyond the study of the present 
paper. In order to compare with other works for the fraction of the 
matter ejected to form the jet in the low luminosity phase (e.g, Yuan 
& Cui 2005), we take the same value of £3 = 0.02 as adopted in 
Yuan & Cui (2005). In a word, although the parameters in the jet are 
uncertain, we still can conclude that the jet is relatively suppressed 
at the high luminosity phase in BHBs. 

The mechanism of jet formation is unclear. Narayan & Mc- 
Clintock (2012) collected a small sample composed of five black 
hole X-ray transients with precise spin measurements. Meanwhile 
the authors estimated the ballistic jet power using the data at 5GHz 
radio observations. It is found that the estimated jet power is cor¬ 
related with the square of a, (with a, = cJ/GM^, J is angular mo¬ 
mentum of the black hole), which is very close to the theoretical 
scaling derived by Blandford & Znajek (1977). However, by sepa¬ 
rately considering the ballistic jet and the hard steady jet, Fender, 
Gallo & Russell (2010) found that there is no evidence for the cor¬ 
relation between the jet power and the black hole spin. Recently, 
an interesting paper argued that the jet formation may be correlated 
with the hot plasma, namely, the jet power is correlated with the 
RIAF when the Fddington ratio is less than ~ 1% and the jet power 
is correlated with the hot corona above the cool disc when the Fd¬ 
dington ratio is greater than ~ 1% (Wu et al. 2013). Our study may 
put some constrains on the mechanism of jet formation, i.e., by sug¬ 
gesting that the relative strength of the jet power may be inversely 
correlated with the Fddington ratio in an accreting black hole (e.g., 
Kbrding, Falcke, & Markoff 2002; Fender, Gallo, & Jonker 2003). 

As a comparison, in the right panel of Figure we also plot 
the relation between Ts-gkcv/f'Edd and ^'S.soHz/i'Edd for GX 339-4 
with the sign of green ‘A’, and V404 Cyg with the sign of or¬ 
ange For Lx > lO^^LEdd, the radio/X-ray correlation with 
f'S.SGHz f'MkcV expected to be interpreted within the frame¬ 
work of the RIAF-jet model, as suggested for the black hole X-ray 
transient XTE Jill 8+480. However, how to justify the existence of 
the RIAF at the high luminosity phase? Qiao & Liu (2009) studied 
the effect of the viscosity parameter a on the critical mass accretion 
rate for the transition from a RIAF to a disc corona system within 
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the framework of the disc evaporation model. The authors derived 
that Merit oc Mndd, i-e., a larger value of a can increase the 
critical mass accretion rate for the existence of the RIAF. A simi¬ 
lar result to the critical mass accretion rate for the existence of the 
RIAF, i.e., Merit cd Mndd was also derived from the RIAF so¬ 
lution (Narayan & Yi 1995b; Mahadevan 1997). By summarizing 
the observational data of dwarf nova outbursts, outbursts of X-ray 
transients, and variability in AGN and so on, King et al. (2007) in¬ 
ferred that the value of a is in the range of ~ 0.1 - 0.4. Extremely, 
Narayan (1996) took a = 1 to explain the BHB systems with high 
transition luminosities. Theoretical understandings for the viscos¬ 
ity have been studied for many years since the pioneering work by 
Shakura & Sunyaev (1973) (Balbus & Hawley 1991; Hawley et al. 
1995; Hawley & Krolik 2001). The estimated value of a from nu¬ 
merical simulations is roughly one order of magnitude smaller than 
that of the observations (Stone 1996; Hirose et al. 2006). However, 
indeed there are numerical simulations showing that the value of a 
sometimes even can exceed unity in the corona (e.g., Machida et al. 
2000 ). 


5 DISCUSSION 

In this work, we proposed a disc corona-jet model to explain the 
observed radio/X-ray correlation of Lr oc for Lx/Tudd <: 10“^ 
in BHBs. We noted that a similar disc corona-jet model was also 
proposed for explaining this radio/X-ray correlation, in which the 
X-ray emission is also dominated by the disc and corona, and the 
radio emission is dominated the jet (Huang et al. 2014). Both in our 
work and Huang et al. (2014), for the disc-corona model, it is as¬ 
sumed that the magnetic field is generated by dynamo action in the 
accretion disc, then due to buoyancy, the magnetic loops emerge 
from the accretion disc into the corona and reconnect with other 
loops. In this way the accretion energy is released in the corona 
as thermal energy and eventually emitted away via inverse Comp¬ 
ton scattering. By studying the energy balance of the disc and the 
corona, Huang et al. (2014) solved the structure of the cold ac¬ 
cretion disc. However, they can not self-consistently determine the 
temperature of the corona T and the Compton scattering optical 
depth T, which are two very important quantities for determining 
the shape of the Compton emergent spectrum. Theoretically, T and 
r should be solved independently for determining the Compton 
emergent spectrum. In Huang et al (2014), it is argued that the ob- 
servationally inferred value of t is in the range of ~ 0.1 - 0 . 8 , then 
through their calculation, t = 0.5 is fixed to solve T for fitting this 
radio/X-ray correlation. In our model, we do not fix the value of 
r. We performed self-consistent Monte carlo simulation to treat the 
structure of the disc-corona system. The fraction of the dissipated 
energy in the corona /, the temperature of the corona T and the 
Compton scattering optical depth t can be self-consistently deter¬ 
mined, which are important advantages of our model. 

We note that in our model, currently in order to simplify 
the calculation of the complex interaction between the disc and 
corona, we always set the albedo, i.e., ‘a=0’, to conduct the Monte 
Carlo simulation, which means that the irradiation photons from 
the corona are fully absorbed by the accretion disc, then are repro¬ 
cessed as the soft photons for the inverse Compton scattering in 
the corona. Since albedo ‘a=0’ is adopted, we don’t have reflec¬ 
tion component in the emergent spectra. We should keep in mind 
that in our model, in the gas pressure dominated case, the origin 
of soft photons for the inverse Compton scattering is dominated 
by the reprocessed soft photons rather than the intrinsic soft pho¬ 


tons of the accretion disc itself. A change the value of the albedo 
‘a’, e.g., an increase the value of ‘a’ means the soft photon lumi¬ 
nosity caused by reprocess decreases, so the Compton luminosity 
decreases, meanwhile the reflection luminosity increases. In our 
Monte Carlo simulation, for the energy conservation, we always 
set that the emergent luminosity is equal to the released gravita¬ 
tional energy. So a change the value of the albedo ‘a’ in a reason¬ 
able range will change the shape of the X-ray spectrum (relative 
strength of the Compton component and the reflection component), 
but will slightly change the X-ray luminosity. So the effects of the 
albedo on the emergent spectrum have only very little change to 
our results. 

As we know, the observed high frequency quasi-periodic os¬ 
cillation (HFQPO), e.g., 150-450 Hz, are consistent with the Kep- 
lerian frequency near the ISCO of a Schwarzschild black hole with 
masses 15-5 solar masses. Meanwhile, the observed HFQPOs do 
not change significantly despite the sizable change in the X-ray lu¬ 
minosity, suggesting the connections between the HFQPOs and the 
mass and spin of the black hole. If the black hole mass is well con¬ 
strained, the HFQPOs can be used to measure the spin of the black 
hole. The observed pairs of frequencies in a 3:2 ratio suggest that 
the HFQPOs are probably produced by some types of resonance 
mechanism (e.g., Abramowicz & Kluzniak 2001). In the present 
disc-corona model, for simplicity, we assume the accretion disc al¬ 
ways extending down to the ISCO of a Schwarzschild black hole, 
i.e, ISCO is fixed at 3Rs- The incorporation of the effects of the 
spin to the disc-corona model is necessary in the future to make the 
model more realistic, meanwhile to match the observed HFQPOs. 

Observationally, there is a positive correlation between the Ed¬ 
dington ratio A and hard X-ray index F for T > 0.01 (Wu &. Gu 
2008, Qiao & Liu 2013). In the present disc-corona model, the 
hard X-ray photon index is F ~ 2.1, and does not change with 
the mass accretion rates. In the present paper, for simplicity, we 
only consider the gas-pressure dominated case for the structure of 
the accretion disc, i.e., only the hard-state solution is considered 
(Liu et al. 2002; 2003). As discussed in Liu et al. (2003), the gas- 
pressure dominated accretion disc can exist for all the mass ac¬ 
cretion rates. When the system is accreting at M > L2MEdd, the 
radiation-pressure dominated accretion disc can exist extending to 
50Rs, which predicts a soft-state solution. When the accretion rate 
is at 0.3MEdd £ M < L2MEdd, the accreting system can be at a 
state between the hard state and the soft state with hard X-ray in¬ 
dex varying with mass accretion rates. 

Meanwhile, in Huang et al. (2014), the jet power is estimated 
according to the proposed hybrid jet model (Meier 1999, 2000), 
from which the radio luminosity is estimated based on the empir¬ 
ical relation from Cyg X-1 and GRS 1915-t-105 (Falcke & Bier- 
mann 1996; Heinz & Sunyeav 2003; Heinz & Grimm 2005). In the 
present paper, based on the internal shock scenario, we calculated 
the emergent spectrum of the jet, which makes our results more 
easy to compare with observations. 


6 CONCLUSION 

In this work, we investigate the radio/X-ray correlation of Lr oc 
for Lx/Undd <: 10 “'^ within the framework of a disc corona-jet 
model, in which a fraction of the matter, rj, is assumed to be ejected 
to form the jet. We calculate the slope of the radio/X-ray correlation 
by assuming a constant 77 for different M. For rj = 0.2%, it is found 
that Lg.sGHz “0 ^‘ 2 ^ 10 -kcw "'ith ^Iq ~ 1.35. For 77 = 0.5%, we derive 
that Lg.sGHz “0 with ^Iq ~ 1.32, which is very close to the 
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Figure 2. 

Left panel: Emergent spectra of the disc corona-jet model around a black hole with M = IOMq assuming a = 0.3 for modeling the radio/X-ray correlation of 
the black hole X-ray transient H1743-322 for ^Edd > 10“^. From the bottom up, the solid lines are the combined emergent spectra of the disc 

corona-jet model for M = 0.02,0.05,0.1,0.3 and 0.5 Mndd: and the corresponding dotted lines are the emergent spectra from the jet with t] = 0.54%, 0.52%, 
0.57%, 0.62% and 0.62% respectively. Right panel: f-s.SGHz/^Edd as a function of Ls-gkcv/^Edd- The red are the observations for H1743-322, the green 
A’are the observations for GX 339-4, and orange are the observations for V404 Cyg. The dotted line is the best-fitting linear regression of H1743-322 for 
Tj-skev/ ^Edd > 10“^. The dashed line is the best-fitting linear regression of GX 339-4 and V404 Cyg. The thick solid line is the model line, and the model 

spectra are shown in the left panel. 


8.5GHz 2 lOkeV 8.5GHZ 2 10KEV 




Figure 3. 


Emergent spectra of the jet for different parameters. In the left panel, M = IOMq, M = 0.5MEdd> V = 0.5%, p = 2.1 and eb = 0.02 are fixed, and from the 
bottom up, Ee are 0.01, 0.04, 0.08 and 0.1 respectively. In the right panel, M = IOMq, M = 0.5MEdd. V = 0.5%, p = 2.1 and Ee = 0.04 are fixed, and from the 

bottom up, eb are 0.01, 0.02, 0.05 and 0.1 respectively. 


case of 77 = 0.2%. As an example, for different M, by changing the 
value of rj, we fit the observed radio/X-ray correlation of black hole 
X-ray transients H1743-322 for is-gkev/^Edd > 10“^. It is found 
that rj is weakly dependent on M, and the mean htting value of 77 is 
~ 0.57%. We note an interesting result, i.e., the mean fitting result 
of 77 0.57% for the radio/X-ray correlation of Lr oc during 

the high luminosity phase is systematically less than that of the case 
for the low luminosity phase (at least 77 > 1 %), which may put some 
constraints on the jet formation, i.e., by suggesting that the strength 


of the jet power is relatively suppressed during the high luminosity 
phase in BHBs. 
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